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Abstract 

Estimating tight resource level bourds is a fundamental 
problem in the construction of flexible plans with resource 
utilization. In this paper we describe an efficient algorithm 
that builds a resource envelope, the tightest possible such 
bound. The algorithm is based on transforming the 
temporal network of resource consuming and producing 
events into a flow network with nodes equal to the events 
and edges equal to the necessary predecessor links between 
events. The incremental solution of a staged maximum 
flow problem on the network is then used to compute the 
time of occurrence and the height of each step of the 
resource envelope profile. The staged algorithm has the 
same computational complexity of solving a maximum flow 
problem on the entire flow network. This makes this 
method computationally feasible for use in the inner loop 
of search-based scheduling algorithms. 

Introduction 

Retaining flexibility in the execution of activity plans is a 
fundamental technique for dealing with the uncertain 
conditions under which the plans will be executed. For 
example, flexible plans allow explici: reasoning about the 
temporal uncontrollability of exogenous events (Morris, 
Muscettola, Vidal 2001) and the incorporation of 
execution countermeasures within the flexible network. 
Tightly constrained schedules (e.g., schedules that assign 
a precise start and end time to all activities) are typically 
brittle and it is very difficult to closely follow their 
directions during execution. For an example of what 
overly tight schedules can do to an intelligent execution 
system, consider the “Skylab strike” (Cooper, 1996), when 
during the Skylab 4 mission astronau:s went on a sit-down 
strike after 45 days of trying to catch up with the demands 
of a fast paced schedule with no room for them to adjust to 
the space environment. 

A major obstacle to building flexible schedules, however, 
remains the difficulty of accurately estimating the amount 
of resources that a flexible plan may need across all of its 
possible executions. This problem is particularly difficult 
for resources with multiple capacity that can be both 
consumed and produced. In the worst case large plans may 
exhibit both a high level of activity parallelism and a large 
number of required synchronizatior constraints among 
activities. Most of the scheduling nethods available to 
date for this problem (Cesta, Oddi Smith, 2000) eventually 
produce a fixed activity schedule, even if they make 
substantial use of an activity plan s flexibility during 
schedule construction. 

To appreciate the difficulty of precisely estimating 
resource consumption, consider the fact that a flexible 
activity plan has an exponential lumber of possible 


instantiated schedules. This means that methods based on 
complete enumeration are typically out of the question. 
Lately, however, new techniques have been developed 
(Laborie, 2001) based on direct propagation of 
information on the temporal constraints of the plan. This 
yields both an upper bound and a lower bound on the 
resource level required by the plan over time. This 
information can be used in various ways, e.g., to decide 
when to backtrack (when the lower/upper bound interval 
is outside of the range of allowed resource levels at some 
time) and when a solution has been achieved (when the 
lower/upper bound interval is inside the range of allowed 
resource levels at all times). Bound tightness is extremely 
important computationally since both as backtracking and 
termination criteria it can save a potentially exponential 
amount of search when compared to a looser bound. 

A natural question is whether constructing the tightest 
possible resource level bounds is computationally feasible. 
This paper answers this question in the affirmative. We 
describe an efficient algorithm for the computation of a 
resource level envelope , a resource level bound such that 
for each time there exists at least a schedule for the 
activity plan that will consume the amount of resource 
indicated by the bound. The algorithm is polynomial, with 
complexity equivalent to solving a maximum flow 
problem on a flow network of the size of the original 
activity plan. 

In the rest of the paper we first introduce the formal model 
of activity networks with resource consumption. Then we 
review the literature on resource contention measures and 
show an example in which the current state of the art in 
resource level bounds is inadequate. Then we give an 
intuitive understanding of our method to compute the 
resource envelope. Then we establish the connection 
between maximum flow problems and finding sets of 
activities that have the optimal contribution to the 
resource envelope. We then show that these sets of 
activities compute an envelope. We then describe an 
efficient envelope algorithm and its complexity. We 
conclude discussing future work. 

Activity Networks and Resource Consumption 

Figure 1 shows an activity network with resource 
allocations. The network has two time variables per 
activity, a start event and an end event (e.g., e ls and e le for 
activity Ai), a non-negative flexible activity duration link 
(e.g., [2, 5] for activity Ai), and flexible separation links 
between events (e.g., [0, 4] from to e 4 S ). A time origin, 
T s , corresponds to time 0 and supports separation links to 
other events. We assume that all events occur after T s and 



before an event T e rigidly connected to T s . The interval 
T s T e is the time horizon T of the ne work. 
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Figure 1: An activity network with resource allocations. 

Time origin, events and links constitute a' Simple 
Temporal Network (STN) (Dechter, Meiri, Pearl 1991). 
Unlike regular STNs, however, each event has an 
associated allocation variable with real domain (e.g., r 3 i 
for event e 3 S ) representing the amount of resource 
allocated when the event occurs We will call this 
augmented network R a piecewise-constant Resource 
allocation STN (cR-STN). In the following we will assume 
that all allocations refer to a single, multi-capacity 
resource. The extension of the results to the case of 
multiple resources is straightforward. An event e* with 
negative allocation is a consume?, while an e + with 
positive allocation is a producer . 

Note that an event can be either a consumer or a producer 
in different instantiations of the allocation variables (e.g., 
event e 2s for which the bound for r?i is [-1, 3]). This 
allows reasoning about dual-use activities (e.g., starting a 
car and running it both make use of the alternator as a 
power consumer or producer). Moreover, some events can 
have opposite resource allocation of other events (e.g., e le 
vs. e ls ). This allows modeling reusab e allocations, such as 
power consumption by an activity. Note that this model 
does not cover continuous accumulation such as change of 
energy stored in a battery over time. A conservative 
approximation can however be achieved by accounting for 
the entire resource usage at the activity start or end. We 
will always assume that the cR-STN is temporally 
consistent. From the STN theory, this means that the 
shor test-path problem associated tc R has a solution. 
Given two events ei and e 2 we denote with |eie 2 | the 
shortest-path from ei to e 2 . We will call a full instantiation 
of the time variables in R a schedule s(.) where s(e) is the 
time of occurrence of event e according to schedule s. We 
will call S the set of all possible consistent schedules for 
R. Each event e has a time bound [ei.(e), lt(e)], with et(e) 

= -|eT s | and lt(e)= |T s e|, representing the range of time 
values s(e) for all se S. Finally, given three events, e 1? e 2 
and e 3 , the triangular inequality |e] e 3 | < |eie 2 | + |e 2 e 3 | 
holds. 

A fundamental data structure used in the rest of the paper 
is the precedence graph, Pred(R), fcr a cR-STN R. This 
is defined as a graph with the same events as R and such 
that for any two events ei and e 2 with |ei e 2 | < 0 there is a 


path from e! to e 2 in Pred(R). Alternatively, we can say 
that an event ei precedes another e 2 in the precedence 
graph if ej cannot be executed before e 2 . There are several 
possible precedence graphs for a network R. A way to 
build one is to run an all-pairs shortest-path algorithm and 
retain only the edges with non-positive shortest distance. 
Smaller graphs can be obtained by eliminating dominated 
edges, e.g., by applying dispatchability minimization 
(Tsamardinos, Muscettola, Morris 1998). ^The cost of 
computing Pred(R) is bound by 0(VE + V 2 lg V) where 
V is the number of events and E the number of temporal 
distance constraints in the original cR-STN. The use of 
different precedence graphs may affect algorithm 
performance but does not affect the theoretical foundation 
described here. 

Considering again the activity network in Figure 1, Figure 
2 depicts one of its precedence graphs with each event 
labeled with the time bound and the maximum allowed 
resource allocation. 

Resource Contention Measures 

Safe execution of a flexible activity networks needs to 
avoid resource contention, i.e., the possibility that for 
some consistent time assignment to the events there is at 
least one time at which the total amount of resource 
allocated is outside the availability bounds. There are 
essentially two methods for estimating resource 
contention: heuristic and exact. Most of the heuristic 
techniques (Sadeh, 1991)(Muscettola, 1994) (Beck et al., 
1997) measure the probability of an activity requesting a 
resource at a certain time. This probability is estimated 
either analytically on a relaxed constraint network or 
stochastically by sampling time assignments on the full 
constraint network. The occurrence probabilities are then 
combined in an aggregate demand on resources over time, 
the contention measure. Probabilistic contention can give 
a measure of likelihood of a conflict occurring. However, 
it is not a safe measure, i.e., the fact that it does not 
identify any conflict does not exclude the possibility that 
the cR-STN could have a variable instantiation with 
inconsistent resource allocation. Exact methods avoid this 
problem and are based on the computation of sufficient 
conditions for the lack of contention. (Laborie, 2001) has a 
good survey of such methods. Current exact methods 
operate on relaxations of the full constraint network. For 
example, edge-finding techniques (Nuijten, 1994) analyze 
how an activity can be scheduled relatively to a subset of 
activities, comparing the sum of all durations with a time 
interval derived from the time bounds of all the activities 
under consideration. Relying only on time bounds ignores 
much of the inter-activity flexible constraints and tend to 
be effective only when the time bounds are relatively tight. 
Therefore algorithms that use these contention measures 
tend to eliminate much of the flexibility in the activity 
network. (Laborie, 2001) goes further in exploiting the 
information about mutual activity constraints. One of the 
two metrics proposed in that paper is the balance 
constraint , an event-centered approach that estimates 
upper and lower bounds on the resource level immediately 
before and after each event e in the cR-STN. These bounds 
precisely estimate the contribution of events that must 



precede e and overestimate the contribution of events that 
may or may not happen before e, The over-estimate is 
obtained considering only the worst-case situation in 
which only the events that have the worst contribution 
(producers for upper bounds and consumers for lower 
bounds) happen before e. Although the balance constraint 
better exploits the information in tf e activity network, the 
bounds that it produce may be very loose for networks 
with significant amounts of parallel sm. 
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Figure 2: Precedence graph with timt /resource usage. 

For example, consider the activity graph in Figure 3, 
consisting of two rigid chains of n activities with the same 
fixed duration and the same fixed activity separation, and 
occurring on a horizon T wide enough to allow any 
feasible ordering among them. Each activity has a 
reusable consumption of one unit and the resource has two 
available units of capacity over time It is clear that all the 
executions of this activity network ^re consistent with the 
resource constraint, since the maximum resource 
consumption is one unit of capacity for each chain at any 
time. However, the balance constraint will always detect 
an over- allocation unless the network is further 

constrained in one of two ways: a) the start activity n of 
one chain occurs no later than the start of the second 
activity of the other; or b) more than two activities overlap 
and there is an activity k on one chain that must start 
between the end of activity i and the start of activity i+2 
on the other chain (Figure 3). These additional constraints 
unnecessarily eliminate a large numl>er of legal executions 
of the activity network. 

The cause of the inability of the balance constraint to 
correctly handle this situation depends on its inability to 
account for the constraint structure of parallel chains 
simultaneously since it can only take advantage of the full 
structure of the network for chains of predecessors. In this 
paper we will show that it is possible to effectively use all 
precedence constraints in the neiwork simultaneously 
leading to the estimate of the best possible upper bound for 
resource consumption. 

Resource Envelopes 

Our approach is to build the tightest possible resource- 
centered exact contention measure. This means that for 
any possible time value we will compute the maximum 
and minimum possible consumption among all possible 
schedules of R. Note that the maximum (minimum) 
overall resource level induced by R for any possible 
schedule can always be obtained by assigning each 
allocation variable to its maximum (minimum) possible 
value. For any specific value assignment to the allocation 
variables, each event has a constant weight: positive, c(e ), 


for a producer and negative, — c(e ), for a consumer. More 
formally, given a schedule seS and a time t € T, E s (t) is 
the set of events e such that s(e) < t. For any subset A of 
the set of events in R, E(R), we will call the resource level 
increment A(A) =2 e +, e -€A c ( e+ ) “ c ( e )• The increment is 
also defined for the empty set as A(0)=O. Therefore, the 
resource level at time t due to schedule s is L s (t) = 
A(E s (t)). The maximum resource envelope at time t is 
L niax (t) = maxses (L s (t)). Similarly, the minimum resource 
envelope at time t is L niin (t) = min seS (L s (t)). Our goal is 
to compute both L imx and L niin over T. 
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Figure 3: Over-constraining of activity network flexibility. 

In the following sections we will rigorously develop an 
efficient algorithm to compute envelopes. Here, we want 
to give an intuitive account of what is involved and of the 
key complexity of envelope computation by analyzing 
some simple examples. 

First consider a single activity with a reusable allocation 
(Figure 4(a)). We could build the envelope L inax by asking 
at each time teT whether Aj can happen before, after or 
can overlap t. If the activity starts with a resource 
production (Figure 4(b)), then we want Ai to start, contain 
or end at t. This is always possible between et(ei s ) and 
lt(e 2s ). During this interval the resource envelope is 1 
while outside of it it is 0. Conversely, if Aj starts with a 
consumer (Figure 4(c)), then we want for the activity to 
happen completely before or after t. This is possible only 
before lt(ei s ) and after et(e le ). The envelope will be 0 at 
every time except between lt(e ls ) and et(eu) where it will 
be -1. This suggests a strategy that looks at each event and 
considers the incremental contribution of the event’s 
weight to the envelope at the earliest time for producers or 
at the latest time for consumers. 

When computing L max for a complex network, however, 
events cannot usually be scheduled independently. 
Consider the simplest network, i.e., a rigidly linked pair of 
activities with a reusable resource allocation (Figure 5(a)). 
In this case the time of occurrence of e 2e and e 3s are bound 
together. Looking at the contribution to the envelope of 
each event in isolation, we would want to add the 
contribution of e 2e as late as possible since it is a 
consumer, and the contribution of e 3s as early as possible, 
since it is a producer. The decision on which time to 
choose depends on the total contribution of both events. 
The total contribution will be added at lt(e 2e ) if tota l 
contribution is a consumption (Figure 5(b)) or at et(e 3s ) if 
the total contribution is a production (Figure 5(c)). Note 
that in both cases e 2s and e 3s are pending at the selected 
time, i.e., their contribution has not been added yet to the 
envelope but they both could occur at the selected time. 
This suggests a strategy that considers all pending events 
at either the earliest time or the latest time of some event 
and schedule those that either must be scheduled or are 
advantageous, i.e., contribute overall with a production of 
resource. 



Now consider the network in Figuie l and the event time 
bounds, maximum resource allocation and precedence 
graph in Figure 2. Assume that we want to compute 
Lmax(3), the maximum envelope at time 3. The set of 
events that may be scheduled before, at or after time 3 is 
{eis, e le , e 3s , e 3e > e 4s }. However, of hese only {ei e , e 3s , e 3e , 
e 4s } are pending since it is advantageous to consider e ls at 
its earliest time 1. The subset of events that we could 
consider at time 3 are all those that will have to occur at or 
before 3 assuming that we select fcr some set of events to 
occur at 3. These subsets are {e 3s }, {e ls }, {e 3e , e 3s } and (e 4s , 
e le , e 3e , e 3s }. Unfortunately, each )f these subsets has a 
negative weight and therefore none of them is considered 
at time 3. At time 4 the set of pending events is 
augmented with e 2s and the total contribution of the new 
subset of pending events {e^, e 4s , ^ e 3e , e 3s} is positive. 
<f0,3],r l > <[5, LO], Ti> 
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Figure 4: Maximum resource envelop.' for a single activity. 

The selection of a maximally advartageous subset among 
the pending events is the key source of complexity of 
envelope calculation. An exhaustive enumeration of all 
subsets can obviously be very expensive. Fortunately we 
can make very good use of the information in the 
precedence graph. It turns out i hat this problem is 
equivalent to a maximum flow p oblem solved on an 
appropriate auxiliary flow network built on the basis of 
Pred(R). We will discuss this rigorously in the rest of the 
paper. 

Calculating Maximum Resource 
Level Increments 

Consider now an interval HcT. We can partition all 
events in R into three sets depencing on their relative 
position with respect to H: 1) the closed events C H with 
all events that must occur strictly before or at the start of 
H, i.e., such that that lt(e) < start(H); 2) the pending 
events R h with all events that can occur within or at the 
end of interval H, i.e., such that lt(ei > start(H) and et(e) 
< end(H); and 3) the open events Oh with all events that 
must occur strictly after H, i.e., such that et(e) > end(H). 
The set R H could contain events that can be scheduled 
both inside and outside H. If H=T, then Ct= 0> Rt = 
E(R) and O H =0- The interval H could be reduced to a 
single instant of time, i.e., H=[t, t] In this case we will 
use the simplifying notation C t =([ t , 1 1 ? Rt=R[t, t] and 

0t=0[ t , tj* 

We are interested in a particular kind of subset of Rh- 
Assume that we wanted to compu e the resource level 
increment for a schedule s at a time teH. This will always 
include the contribution of all events in Ch and none of 


those in O h irrespective of s and t. With respect to the 
events in R H , we can see that if an event is scheduled to 
occur at or before t then all of its predecessors (according 
to Pred(R)) will also have to occur at or before t. In other 
words, it is possible to find a set of events X e Rh such 
that the events e p eRH that are scheduled no later than t in 
s are those such that |e x e p | < 0 for some e x e X. We call 
this the predecessor set of X, Px« Therefore, the resource 
level at time t for a given schedule s is the sum of the 
weights of events in Ch and in Px- 
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Figure 5: Maximum envelope for two chained activities. 

It is easy to verify that given two predecessor sets Px and 
P Y , both P x n Py and P x uP Y are also predecessor sets. 

Resource Level Increments and Maximum Flow 

Since we are interested in the maximum resource level we 
want to find the predecessor set with maximum resource 
level increment. We will do so by finding a maximum 
flow for an auxiliary flow network built from Rh and 

Pred(R). 

Resource Increment Flow Problem: Given a set of 
pending events Rh for cl cR-STN R y we define the resource 
increment flow problem F(Rh) with source (7 and sink tcls 
follows : 

1. For each event e e R H there is a corresponding node 

eEFfRn)- + 

2. For each event e+ eRh, there is an edge e with 
capacity c(e + ). 

3. For each event e e Rh, there is an edge e'->T with 
capacity c(e), i.e., the opposite of e~ 's weight in R . 

4. For each pair ofej and ej with an edge e\—>e 2 in the 
precedence graph Pred(R), there is a corresponding 
link e l —*> 2 in F(R U ) with capacity +<~. 
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Figure 6: Resource increment flow problem. 



As an example. Figure 6 shows the auxiliary flow problem 
for Rt relative to the activity network in Figure 1. 

A detailed discussion of flow problems is beyond the scope 
of this paper (for a complete treatment see (Cormen, 
Leiserson, Rivest 1990). Here we highlight some 
fundamental concepts and relations that we will use. We 
will indicate as f(ei, e 2 ) the flow associated to a link 
ei _->e 2 in F(R h ). The flow function is skew-symmetric, 
i.e., f(e 2 , e 2 ) = - f(e b e 2 ). Each flow has to be not greater 
than the capacity of the link to which it is associated. For 
example, referring to the flow network in Figure 6, 0 < f(a, 
e 2e ) < 2, 0 < f(ei e , x) < 4 and f(e 4e e 4s ) > 0. Note that a 
flow from ei to e 2 can be negative only if the flow network 
contains an edge e 2 — >ei with posit. ve capacity. We also 
use an implicit summation notation f(A, B), where B and 
A are disjoint event sets in F(Rh)> to indicate the flow f(A, 
B) = SaeA^beB^a, b). Consider now any subset of events 
AcRh and let us call Atfie set of events A = Ru-A. The 
following flow balance constraint always holds: f({o}, A) 
= f(A, {t} ) + f(A, A). The total network flow is defined as 
f({a}, Rh) = f(Rm {x}). The maximum flow of a network 
is a function f inax such that the total network flow is 
maximum. 

The fundamental concept used by all known maximum 
flow algorithms is the residual network . This is a flow 
network with an edge for each pair )f nodes in F(Rh) for 
which the residual capacity, i.e., the difference between 
edge capacity and flow, is positive. Each edge in the 
residual network has capacity equal to the residual 
capacity. For example, considering he network in Figure 
6, assume that f(ei e , x) = 3 and f (a, e 2e ) = 2. The residua) 
network for that flow will have the flowing edges: €i c — >x 
with capacity 1, x— >e ie with capacity 3 and e 2e -^o with 
capacity 2. Also note that any residual network for any 
flow of F(Rh) will always have an edge of infinite capacity 
for each edge in the precedence graph Pred(R). 

In this paper we will make use of three different kinds of 
paths in the residual network. The f rst is an augmenting 
path connecting a to x. The existence of an augmenting 
path indicates that additional flow c m be pushed from a 
to x. Several maximum flow algorithms operate by 
searching for augmenting paths. Alternatively, the lack of 
an augmenting path is the condition that indicates that a 
flow is a maximum flow. The second kind of path is a 
flow-shifting path , a loop connecting r to x which does not 
affect the overall flow in the network. Finally, the third 
kind is a reducing path , i.e., a path from x to a. Pushing 
flow through a reducing path reduces a network’s flow. 

We now establish the relation between the resource level 
increment A(A) and any flow in F'Rh). We define the 
producer weight in A as c(A + ) = 2! e + e a c(e + ) and the 
consumer weight in A as c(A’) = X e - e a c ( e )• a ^ so 
define the producer residual in A as r(A*) = c(A + ) - 
f({o}, A), i.e., the total residual capacity of the edge 
incoming A from s, and the consumer residual in A as 
r(A‘) = c(A') - f(A, {x}). The followir g relation holds. 

Lemma 1: A(A) = r(A + ) - r(A') + f(A, A). 

Proof: A(A) = c(A + ) - c(A‘) = (c(A + ) - f({a}, A)) - (c (A' 

) - f({o}, A)) = r(A + ) - (c(A*) - f(A, {x}) - f(A, A)) = 
r(A + ) - r(A’) + f(A, A). □ 


We now focus on predecessor sets such as P x . 

Lemma 2: f(P x , Lx) ^ Moreover, f(P X) Lx)=0 if and 
only iff(e 2 , e 2 )=0 for each e^Lx and e 2 eP x . 

Proof: From the definition of predecessor there is no edge 
e 2 _> ei in F(R h ) with ejGPx and e 2 eP x Therefore, f(e 2 , 
e {) < o and f(P x , P x) ^ 0* The second condition can be 
demonstrated by observing that the sum of any number of 
non-positive numbers is 0 if and only if each number is 0. 
□ 

Another way to express Lemma 2 is that f(P Xj Ex) = 0 if 
and only if there is no link ei->e 2 in the residual network 
where e^Px and e 2 e P x . 

Corollary 1: A(P X ) <r(P x v ) - r(P x '). 

Proof: Immediate from Lemma 1 and Lemma 2. 

Maximum flows and maximum resource level 
increments 

We are now ready to find the maximum positive resource 
level increment. Note that we are not interested in event 
sets with negative resource increments since, as we 
discussed before, we will only account for events in our 
resource envelope simulation if they have a positive 
contribution. If they do not, we will take them into account 
when we must, i.e., when their temporal upper bound 
becomes lower or equal to the current simulation time. 

First we address the problem of whether Rh contains a set 
of predecessors P* with positive resource level increment, 
i.e., A(P*) > 0. To do so we will make use of a maximum 
flow f nulx of F(Rh). We will indicate with r inax (A) 
producer/consumer residual computed for f rn ax- The 
following fundamental theorem holds. 

Theorem 1: Given a partial plan Rh, there is a 
predecessor set P* such that A(P*)>0 if and only if 
r,nax(PH + ) > 0 . 

Proof: D : We prove that if there is a P such that A(P )> 

0. then r imx (R H + ) > 0. Assume that r inax (R H *) = 0. This 
means that for any predecessors subset P X we have 
r uiax (Px") = 0. From Lemma 2 we would have A(P ) < 
-r inax (P*') < 0, that is a contradiction. 

We prove that if r nwx (R H + ) > 0 then we can identify a 
P* such that A(P*)>0. If r niax (R H + ) is positive, there must 
be some e + such that r niax (e + ) > 0. Let’s select P as the set 
of events reachable by some path in the residual network 
originating from e + . The following three properties hold. 

1 . P * is a predecessor set. 

If not, there will be an event e 2 £ P such that |eie 2 | < 0 
for some event ei£P . From the definition of Pred(R), 
however, we know that there must be a path in Pred(R) 
from ei to e 2 . Since this path will be present in F(R H ) 
with all links having infinite capacity, the path will also 
always be present in any residual network for any flow. 
Therefore there is a path in the residual network going 
from e + to e 2 to e 2 and e 2 €P", which is a contradiction. 

2. r max (P*‘) = 0. 

If not, there will be an event e e P such that r lliax (e ) > 
0. We can therefore build an augmenting path of F(R H ) 
as follows: 1) an edge a->e + with positive residual 
capacity r inax (e + ); 2) a path in the residual network from 



e + to e", which exists by definition of P*; and 3) an edge 
e'->x with residual capacity iw(e'). The existence of 
the augmenting path means that f niax is not a maximum 
flow, which is a contradiction. 

3.f max (P\ P*) = 0 

Since P* is a predecessor set, from the proof of Lemma 
2 we know that f ma *(P\ El) <0 If fmax(P*> El)<0, it 
means that there is a pair of events eiG P and e 2 gP 
such that f nia .x(ei, e 2 ) < 0. This means that the residual 
capacity from ej to e 2 is positive and therefore there is 
an edge ei— »e 2 in the residual network. But this means 
that e 2 eP\ which is a contradiction. 

Applying the properties of P* to the reflation in Lemma 1 
we obtain A(P*) = r^ (P + ) - r nux (P ) + fma*(P ? E_) = 
r™x(P* + ) ^ r niax (e + ) > 0. □ 

It is now easy to find the predecessor set P,»a X with 
maximum positive resource level increment. 

Theorem 2: Consider all event . e + t eRh such that 
rmax( e +i) > 0 aR d consider the event set P max t-4+i P i 
where P\ is the set of events reac hable from e* in the 
residual network of f^. Then P ma x is a set of 
predecessors in Rh with maximum L\(P max ) > 0- 
Proof: Each of the P*i has the properties proved in 
Theorem 1. We show that P™* also has those properties. 

1. Being the union of predecessor sets, Pmax is also a 

predecessor set. ^ 

2. We know that r niax (e') = 0 for each e £ P l Therefore 

TnaxC-P max) = * * 

3. From Lemma 2 we know that since f(P P j)=0, there is 
no flow from events in P*i to e\ents in P Therefore 
there is also no flow from events in P niax = r\ Ej to 
events in P n ^ x - Hence, from Lemma 2 f(P ma X > 

E*™)=0. 

Therefore A(Pniax) = r llia x(P + max) >0. 

Moreover, since by construction P miX contains all e + i with 
r max (e + i) > 0. Therefore, for any other predecessor set Px it 
is r niax (P + x ) ^ r luax (P + ]nax ). Hence, A(Px) ^ r uiax (P x ) — 
r n iax(P x) — r max(P X ) — **max(P max) = A(Pmax)« 0 

So far we have constructed P nwx from a specific maximum 
flow for F(Rh). However, it turns cut that Pmax is unique 
for all maximum flows of F(R H ). Moreover P ]liax contains 
the minimum number of events among all predecessor sets 
with maximum positive resource lev *1 increment. 

Theorem 3: For any solution of the maximum flow 
problem for F(R H ), P ma x is the minimal predecessor set 
with maximum resource level increment A(P max f 
Proof: Consider the set {f niax4 } with i=l, n of all n 
different maximum flows of F(R H ). Since each 

corresponding P nmX|i is a maximum positive resource level 
increment sets, A(P niax ,i) = A, mx . Also, given two distinct 
maximum flows i and k, we have Pmax,j - Pink u Pk-i 
where P ink = P max ,i n Pmax,k and P k i ^ P™** - Pmax,i- In the 
following we will indicate with rj(e) the residual r niax (e) 
computed in flow f niaX j. 

First we observe that for any two dis inct i and k, A(P in k)= 
A max* In fact, P ink is a predecessor set and rj(P ink ’)= 
r k (Pink')=^- Therefore A(P ink )= ri(Pi nk + )= r k (Pink )• We 
need now to show that rj(P iok + ) = r (Pmax,i + ) - ri(P, mx . k ) • 


In fact, it must be rrfPuX) since P max ,i must contain all 
e + events with rj(e + ) > 0. Also, A niax = n(Pmax,i )=A(P, mx , k ) 
< « < f i (P lnaJ[t k + ) which implies 

ri(Pi-k + )=0. But this means that the flow in A(P iok ) is selt 
contained, i.e., there is no edge in the residual network of 
flow f nwx ,i that exits P ink . Therefore, in this flow none of 
the events in P ik is reachable from an e + event and 
therefore Pi_ k =0. With a symmetric argument we can see 
that Pu-i=0* Therefore for any i and k it must be Pma X ,i= 
P lliax k = P lliax The minimality of P^x derives from 
applying to P inax - P 0 the same argument used to 
demonstrate that Pi. k is empty, where P cPma* is a 
predecessor graph with maximum positive resource level 
increment. □ 

Building Resource Envelopes 

So far we know that the resource level at time t e H for a 
given schedule s is the sum of the weights of the events in 
C H plus those of the events in some predecessor set P x . It 
is not immediately obvious that the converse also applies, 
i.e., that given any predecessor set Px one can determine a 
time t x e H, the separation time , and a schedule s x , the 
separation schedule , such that all and only the events in 
C H uP x are scheduled at or before time t x . The reason this 
is not obvious is that events are still constrained by upper 
bound constraints, i.e., the metric links that are not 
included in Pred(R). Scheduling some event too early 
with respect to t x may therefore force some event to occur 
before time t x whether the event is a successor in Pred(R) 
or not. We will show that indeed we can find a separation 
time and schedule for any P x and therefore also for Pmax- 
For the latter we will show that t x represents one of the 
times at which the resource level is maximum over H for 
any schedule. This will yield the resource envelope L niax it 
we reduce H to a time t and scan t over the horizon T. 

Latest events 

The first step is to identify the events in P x that will be 
scheduled at time t x . We say that e is a latest event of P x 
if it is not a strict predecessor of any other event in P x , 
i.e., for any e P x , |e x e| > 0. There must be at least one 
latest event in P x . If not, for every event e k e P x , there 
would be an event ej e P x such that |e s e k | < 0. But this 
would mean that it would be possible to create a cycle ej(D 
-> e i(2 ) -> ... -> e i(n ) linked by links |e i(k) e i(k+ i)| < 0, which 
is a contradiction to the hypothesis of temporal 
consistency of R. We will call P x , late the set of all latest 
events in P x . Also, we define P x, early = Px - Px,iatc* 

The following properties hold for the temporal relations 
between events in P x , i a te> Px, early an d Ex- 

Property 1: For any two events e h e 2 e Px, iate> l e i e 2 \ 
and \e 2 >0. 

Property 2: For any two events ej € Px an d e 2 Px, iate> 

| e 2 e i\ > 0. 

If not, ei would belong to P x by the definition of P x - 

Property 3: Any event ejE Px, early is o. strict predecessor 
of some e 2 E Px,i<ue> ce., \e 2 e 2 \ < 0. 



If not, consider any two events e 1? e 3 e Px, early For any e 2 it 
would be |e 2 ei|=0 and |e2e 3 |=0. Therefore, 0=|e 2 e3| < \e^ei\ 
+ |eie 3 | = |e 1 e 3 |, i.e., [e^l > 0. Since this would be true for 
any pair of events in Px, early and ter all distances between 
any event in P x , e nriy and any event in Pxjatc, all events in 
Pxjearly WOUld be latest events, i.e., F X ,early=0- 

Separation Time for Latest Events 

We now show how to construct the separation time t x at 
which we will schedule all latest events. 

Lemma 3: There is a time inter\>al t X ,mad that is in 

common among all time bounds iet(e) f lt(e)] with e e 
Pxjate an( l such that start(H) <t x max 
Proof: First, we show that there must be a time value in 
common among all time bounds, f not, there would be 
two events e x , e 2 e P x ,iate such that et(e x ) > lt(e 2 ). From 
the triangular inequality we also have that le^l < - et(e x ) 
+ lt(e 2 ) < 0, which is inconsistent with Property 1. Now, 
assume start(H) > t x , max* By the way the interval [t x , mim 
tx, max] is constructed, there must exist an event e e Palate 
such that lt(e) = t x , mn. For this event it would be lt(e) < 
start(H) that is a contradiction with e belonging to Rh*E 

We define t x = max (start(H), tx,mmh with t x = start(H) if 
p x =0, in which case t x = start (H). We can then show 
that the time bound of each event in P x indicates that each 
of them can be scheduled after t x . 

Lemma 4: For any event e e d(''i) > h 
Proof: By definition of R H it must be lt(e) > start(H). So 
we only need to consider the case in which t x = t X) nun > 
start(H), In this case there is at leas: one event ei e Palate 
such that et(ei) = t x ,mm* For this event it is |ei e| < - et(e t ) 
+ lt(e). From Property 2 we krow that |ei e| > 0. 
Therefore, lt(e) > et(e x ) + |e x e| > et(e x ) = t x ,mi «.0 

Separation schedule for predecessors 

We now show- how to build a separation schedule s x for P x 
and t x , i.e., a schedule such that s x <e) ^ tx for eeCHuPx 
and s x (e) > t x for eePxuO x . Ncte that the following 
discussion holds also if P x =0. 

We will do this with the following algorithm. 

1. Schedule all e’e P x , j a te at t x , i.e., s x (0 = t x . 

2. Propagate time through R obtai iing new time bounds 
[et’(e), lt’(e)] for each eeE(R). 

3. Schedule all events e”e E(R) -P x ,iatc at their new 
latest time, i.e., s x (e”) = It* (e”)- 

To show that s x is a schedule we need to see that it is 
consistent with respect to R. We see that step 1 is 
consistent since: 1) t x belongs to I he intersection of all 
latest event time bounds; 2) since for any pair of latest 
events |eie 2 |>0, scheduling one at t x does not prevent any 
other latest events to be scheduled a time t x as well. Step 
3 above is also consistent since it is always possible to 
schedule all events at their latest times without temporal 
repropagation. 

Now we need to show that the property defining a 
separation schedule is satisfied for s x . Note that we 
already know that it is satisfied for events in P x ,iate- By 
definition of Ch and Oh? we also know that it is satisfied 


for events in these two sets. Therefore, we need to show 
that it is satisfied for P x ,eariy and P x . 

a) It’(e) < t x for all e e Px, ca n y 

According to Property 3 we can pick an event eieP x ,i a(e 
that |e L e| < 0. From the triangular inequality we have 
lt’(e) < It’(ei) + |e x e| < lt’(eO = t x . 

b) lt’(e) > tx for all ee Px- 

From Lemma 4 we know that before the re-propagation 
it was lt(e) > t x . After the propagation, either 
It’(e>=lt(e), in which case the condition is satisfied, or 
It’ (el has changed with a propagation starting from 
some event e x e P x , So it must be lb’(e) = t x + |eie| 
and since from Property 2 |e x e| > 0, lt’(e) > t x . 

We can now compute L luflX over the entire time horizon T. 
Pi«»x(Rh) indicates that it is computed over F(R H ). 

Theorem 4: The maximum resource consumption for any 
schedule of R over an interval HcfT is given by A (C H ) + 

A(P max (R H ))- 

Proof: We know that at any time teH the events in K H 
that are scheduled before t are a predecessor set P x . For 
the resource level at time t it is always A (C H ) + A(Px) < A 
(C H ) + A(Pmax(Rn) ) » latter being the resource level at 
the separation time for the separation schedule. This is 
true also if P™x(Rh) is empty.T 

There are two interesting special cases of Theorem 4. 

Corollary 2: The maximum possible resource 

consumption for R over T is equal to A(P max (Rr)f 

This means that estimating the maximum possible 
resource consumption over time has the same complexity 
of a maximum flow problem. 

Corollary 3: L max (t) = A(C r )+A(Pmax(Rt))‘ 

The last formula tells us how to compute the resource 
envelope. We now' need to find an efficient algorithm. 

Efficient Computation of Resource Envelopes 

From Corollary 3, the naive approach to compute a 
resource envelope would be to iterate over all possible 
teT. We can improve the running time by considering 
that we only need to compute L imx at times when either C t 
or Rt changes. It is easy to see that this can only happen at 
et(e) or lt(e) for any eeE(R). Therefore we need to re- 
compute a maximum flow for a partial network in R only 
2n times, a substantial improvement over |T|. 

The complexity of some known flow algorithms is 
described in Table ](Cormen, Leiserson, Rivest 1990). 
Note that the number of edges E is 0(V*) where V is the 
number of events and l<x<2. Therefore the^ complexity of 
maximum flow algorithm is always 0(V lg*(V)) with 

l<k<5 and j g{ 0, 1}. Let us now consider the worst case 
complexity of re-computing a flow at et(e)/lt(e). In the 
worst case, C t will remain empty at all times and the size 
of each R t will increase by 1 for each computation of the 
flow Therefore the worst case complexity of this method 

is 0(I i= i v O(i k lg i (i))=0(V k+1 lg i (V))? a polynomial of 

higher order than maximum flow. 



We can do better. Assume sorting all earliest and latest 
times in ascending order to yield a set {t(l), t(2), 
t(2n)}. Suppose now that when we compute the maximum 
flow for F(R t( i)) we make as much use as possible of the 
maximum flow for F(R«h))- In this case we can come up 
with an algorithm with the same worst-case complexity as 
computing the maximum flow on th^ entire network. 

Incremental Change of Pending Events 

Before we introduce the algorithm, let us consider the 
differences between R t <i) and R t( M) (l<i<n). The first 
difference 8(C t(i) ) = R t( i-i) - R t a> is the sets of events e such 
that t(i) = lt(e). They must move from R t( i-i> to Ctro at time 
t(i). The second difference, 5(R l(i) ) = R t (o - Rta-n are the 
events e such that t(i) = et(e). They must move from (Vu 
to R t (i) at time t(i). 

Figure 7 gives a complete picture of how all relevant event 
sets change at time t(i). In this picture E inax (t(i-1)) = Ckm) 
u P ITiax (Rt(i-i)) is the set of events reeded to compute the 
resource envelope at time t(i-l). E nwx (t(i)) = C t(i) u 
Pmax(Rt(i)) is used to compute the resource envelope at time 
t(i), with C t (i) = Ctfi-i) + 8(Ct(i)). The difference E inax (t(i)) 
- E nlflx (t(i-l)) can be separated into two disjoint sets 
5(C t {i))- Pmax(Pt(i-l)) and Pniax(Rt{i)) “ Pmax(Rt(i-l))* The goal 
of the efficient envelope algorithm is to identify the set 
Pnax(Rtti)) - Pn»x(Rt(H)) with less effort than computing 
PmarCRtci)) and P ma x(Rt(i-i)) with separate maximum flow 
computation and then differentiating them. 


Ford-Fulkerson 

0(E f\ ) 

Edmonds-Kaip 

O(VF') 

Simple preflow-push 

0(V J E) 

Preflow-push 

0 (Vh 

Goldberg-Tarjan 

0(Vf: lgCV 2 /E)) 


Table I: Complexity of known maximum flow algorithms 

Before proceeding notice that Figure 7 assumes that 
Pmax(Rtti)) n Pn»x(Rt(M,) C P™x(Rta;). As we will see this 
is indeed the case. The consequence of this is that as soon 
as Pnwx(Rtfi-i)) has been determined and accounted for in 
the envelope calculation, the subnetwork of F(R t {i.u) 
consisting of all events in Pmax(R«i n) and all incoming 
and outgoing edges can be deleted. This allows the 
computation of P,mx(Rt<i))-*Pniax(Rt<i-i)) directly from the 
maximum flow of F(R t a)-P ntax which can save 

significant work. 

The second efficiency improvement is computing the 
maximum flow of F(Rt(j)-Pniax(Rt<i-u) by incrementally 
modifying the flow of F(R t {M)'Pnifx(Rt(i- 2 ))) during the 
deletion of events in 8(C t( i)) and the addition of events in 
8(R t( i)) while maintaining the maximum flow property. 

Incremental Modification of Maximum Flow 

Let us focus on modifications of tf e flow network that 
preserve the maximum-flow. To do so we introduce the 
concept of a prefix and postfix of ,i resource increment 
flow network F. Consider a partition of events in the 
network in two event sets, Post(F) and Pref(F). We say 
that Post(F) is a postfix of F and PrtfiF) is a postfix of F 
if for each pair of events eiePost(F) and e 2 EPref(F), |©2 


Ci | > 0. It is immediate to see that for any flow of F it can 
only be f(e 2 , ei)<0. Therefore the residual network 
contains an edge e 2 -+ei only if there is an edge ei— >e 2 in 
the flow network and there is a positive fie!, e 2 ) passing 
through it. 

We can see that 8(C t(i) ) is a prefix of F(A) where A is any 
subset of Rt(M) that contains SCCkd). In fact, consider a 
pair of events e 2 e 8(Ct(i)) and e 2 £ A - 8(C t (i))- From the 
definition of S(C t (,)) we have lt(e 2 ) = t and It(ei) > t+1. 
From the triangular inequality lt(ej < lt(e 2 ) + |e 2 we 
can deduce |e 2 ej > lt(ei) - lt(e 2 ) > t +1- t = 1 > 0. A 
similar argument applies to demonstrate that S(Rtu)) is a 
postfix of F(B) where B is any subset Rt(i) that contains 
5(R t (i)F 


We now introduce two flow modification operations: flow 
reduction and flow expansion . 


Ctfj.i j 


Rlfi-lj 
Rlfi) 

6(C Ul) ) 

5(R, m ) 

P«*(R 
Pmn, ( Rk,)) 

E^CtCi-D) 

Emoiftfi)) 

Figure 7 : Incremental change for set of pending events 

Flow contraction: Consider a network F(A), a flow f for 
F(A) and a prefix of A, Pref(A). Flow contraction consists 
of the following two steps: 

])while there is a flow-shifting path in the residual 
network connecting an e' event in PrefiA) to an e 
event in A-Pref(A) y push flow along the path and 
update the residual network accordingly; 

2)while there is a reducing path in the residual network 
connecting an e event m PrefiA ) to an e event in A- 
PrefiA ) y push flow along the path and update the 
residual network accordingly. 

Lemma 5: If the flow f max is maximum for F(A), flow 
contraction produces a maximum flow for F(A-PrefiA)). 
Proof: Jf fmax is maximum, the flow f produced at the end 
of step L is still a maximum flow for F(A). This because 
flow shifting does not affect f({s), A + ) since no flow is 
pushed back through any edge f(s, e + ). At the end of step 2 
we will have a flow f\ Note however that any 
modification of the residual network in step 2 can only 
eliminate existing edges and therefore eliminate paths. 
Since f maximally implies that there are no augmenting 
paths in it going to events e e A-PrefiA), there will be no 
such augmenting paths in f” either. Therefore f’ will be 
maximum in F(A- PrefiA)). □ 

Note that in achieving the maximum flow for F(A- 
PrefiA)) it is always better to us flow-shifting paths before 
reducing paths, since if flow needs to be moved from 





Pref(A) to A-Pref(A) to achieve 1 optimality, a flow- 
shifting path is always shorter than the concatenation of a 
reducing and an augmenting path. 

Flow expansion : Consider a network F(A), a postfix of 
A } Post(A)), a flow f for F(A-Posi(A)f Flow expansion 
consists of the following steps: 

• while there is an augmenting path in the residual 
network of F( A) connecting an e * event in Post (A) to 
an e event in A, push flow along the path and update 
the residual network accordingly. 

It is clear that flow expansion prod ices a maximum flow 
finax for F(A). Also, if the starting flew for F(A-Pref(A)) is 
a maximum flow, flow augmentation will minimize work. 
This is important in our application of maximum flow to a 
sequence of n F(R lfkJ ) problems since re-doing 
unnecessary work may negatively impact asymptotic 
complexity. 

Flow Network Reduction 

Now we will show that it is always safe to eliminate any 
events in P inax (Rta-i)) from considera ion once their impact 
on Lmax(t(i-1)) has been recorded. Consider the set 
5(C t (i))nP maK (Rt(i.i>). The effect of these events is already 
included in L max (t(i-1)) and therefore their contribution 
does not need to be included in L iu;ix (t(i)). Let us now 
consider the effect that flow reduc ion applied to these 
events has on the maximum flow of F(R t (i_i) n R t (i>)» From 
the property of the predecessor set with maximum 
resource level increment, we know that f(ei,e 2 )=0 and 
f(e 2? ei) = 0 with e^Pn^CRto-n) and £ 2 ^Emax(Rt(i-i))« 
Therefore, the residual network for the maximum flow of 
F(R«i-i)) does not have any edges e 2 ->ei. Also, since all e 
g PniaxQRt(i-u) are saturated (r 1Bax (e')=0), there cannot be 
any edge e 2 — >t in the residual network. Therefore, there 
cannot be any flow shifting and the only way to push flow 
back ffom an e'G P n ux(Rt(i-i)X is to co so through an e + G 
Pniax(Rtft-i))* Therefore, flow changes during flow 
contraction due to events in P m ax(Rtu-i)) do not affect the 
rest of the network. 

Consider now P nW x(Rt{i-i))-5(C tf i)). We know that after 
contraction the flow is maximum for F(R t(i .i)nRt(i)). 
Therefore all e'G P imx (R t ti-i))-8( are still saturated. 
Also, there will be at least one e + e P,nax(Rtfi-i))“5(C t (i)) 
with r niax (e + )>0. If we now add 8(R t <i)) and apply flow 
expansion, we know that throughout the process there will 
be no augmenting path exiting ffom e*G P nwx (R t (i. 
D)-5(C t (i)). At the end of the process, P n *x(Rt(i-i))-5(C t a)) 
will still be isolated and therefore A(P imx (R t(i . 
i ) )-5(C t(i )))>0. Therefore P niax (Rifi))=?Pniax(Rt(i-i))~5(Cui)). 
Therefore, A(P niax (R t a-i))) will be pari of both L nwx (t(i-1)) 
and L llulx (t(i)). Moreover, the presence of P^Rio-d) has 
no effect on the flow modification operations and therefore 
all events in P^Rta-i)) together with their incoming and 
outgoing edges can be safely eliminated from any flow 
network considered by the algorithm after F(R t(i _i)). 

Figure 8 shows the pseudocode of the algorithm. Here we 
assume that after executing flow contraction , the function 
Flow_Contraction function also deletes from the network F 
the portion pertaining to the events in p.Eciose. Similarly 
we expect FIow_Expansion to add to F the portion 


pertaining to Epend. The function Network_Reduction 
deletes from the network F the portion pertaining to Pmax. 
Finally, Extract_P_max finds the P nW x in the maximum 
flow of F by collecting the events that are reachable in the 
residual network of F ffom each e + with r imx (e + )>0. 

Complexity Analysis 

Now we can show that the complexity of 
Resource_Envelope is asymptotically the same as running 
a flow algorithm over the flow network F(R) for the entire 
cR-STN. We demonstrate this for the Ford-Fulkerson 
method. This is the simplest maximum flow method and 
tends not to perform well when the internal “pipes” of the 
flow network have bottlenecks (Cormen, Leiserson, Rivest 
1990). However, this is not true for F(R) in which the 
capacity of all internal pipes is +°°. Therefore in our case 
the Ford-Fulkerson method should perform well. The 
complexity of Ford-Fulkerson is 0(E |f inax |), where E is 
the number of edges in the flow network and |fmax| is the 
amount of the maximum flow pushed through the 
network. 

Resource_Envelope (R, Pred(R)) 

{ Laid := 0; I* envelope level at previous iteration time 7 

Pm«« := 0 /* temporary variable to hold Pm»* V 

Envelope =0; I* envelope profile, a list of pairs e=<t, L> where t is the 
time at which the envelope reaches level L. The 
envelope stays constant at L until the time of the 
following entry in the list 7 

F := 0; f* auxiliary flow graph for F(Rt<o-Pm«(Rt(i-n)V 

E := (sorted list in ascending order of t of triples <Edo«, E p . n d, t> where E C io»« 
is the set of events e in Pred(R) such that lt(e) - t and E P *nd is the list 
of events e in Pred(R) such that et(e) = t}; 

t := 0 

while (E is not empty) 

{ E C ur := 0; 

P := Pop(E); 

t := p.t; 

F := Flow_Contraction{F, p.Eciow); 

F := Flow_Expansion(F, p.E p .nd, Pred(R)); 

Pmi* := Extract_P_max (F); 

F := Network_Reduction(F, P max) I 

Lfi«w 1= Laid + A(Eclou) + A(P max); 

Envelope := append (<t, L n .w>, Envelope); 

Lold Lnawl 

} 

} 

Figure 8: Maximum resource envelope algorithm . 

Let us consider separately the total cost of the three flow 
propagations in Flow_Expansion and Flow_Contraction: the 
one across augmenting paths, the one across flow-shifting 
paths and the one across reducing paths. 

Consider the propagation across augmenting paths. 
Consider the total flow that we will be able to push ffom 
the source throughout these propagations. Since when we 
push flow from e + events in 5(Rm)), some events e' may 
have become unavailable (through deletion) or saturated 
(due to flow shifting). Therefore the total flow we can 
push in Fiow_Contraction is no greater than |f niax | in F(R). 
Moreover, the set of edges over which the search for 
augmenting paths is conducted is always not greater in 



each invocation of Flow_Expansion than in the the 
application of maximum flow to the entire F(R). 
Therefore, the total cost of Flow_Expan$ion is 0(E |f nwx |) 
Consider now the propagation across flow-shifting paths 
in F!ow_Contraction. Consider an auxiliary flow graph 
Shift(R) built as follow: for each 5(C t(i )) consider all edges 
in the flow-shifting paths before the application of 
Fiow_Contraction(F, 5(Ci (i ))) and for each pair of events e x 
and e 2 (including ct and x) add the link with maximum 
capacity in some of the flow-shitting paths across all 
invocations of Flow_Contraction. The total number of 
edges in Shift(R) is bounded by 2 E + 2 V, since this is 
the maximum number of edges n a residual network. 
Considering now pushing flow from an e' in 5(C U i)), we 
notice that the pipes that are available in Shift(R) are no 
smaller than those available in the flow-shifting paths 
during the execution of Flow_Contraction. Moreover, the 
lengths of the paths is no shorter in Shift(R) than in each 
invocation of Flow_Contraction since events are not deleted 
in Shift(R). Finally, the total flow that can be shifted is no 
greater than |f IIMX |, since this is an upper bound of the 
maximum amount of flow that reaches x from some e . 
From this argument we can see that the total amount of 
work done during the flow-shifting step F!ow_Contraction 
is bounded by 0(E |f inax |). A similar argument applies to 
the flow reduction step in Flow, Contraction with its 
complexity again bound by 0(E |f n iax|)* Therefore the total 
cost of Flow_Contraction and Flow_Expansion in 
Resource_Envelope is bounded by 3 0(E |finax|)- 
Considering the other steps in Resource.Envelope, the 
sorting step to initialize E is 0(V IgV), and the total cost 
of Extract_P_max and of lncremenlally constructing and 
deleting the flow network, is 3 0(E). Under the reasonable 
assumption that the asymptotic cosi of flow computation 
dominates 0(V IgV), the total cost of Resource_Envelope 
is 0(E |f inax |), i.e., it has the same asymptotic complexity 
than running the flow algorithm once over the entire 
F(R). 
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Conclusions 

In this paper we describe an efficient algorithm to 
compute the tightest exact bound on the resource level 
induced by a flexible activity plan. This can potentially 
save exponential amounts of work with respect to looser 
bound computations. Future work includes testing the 
practical effectiveness of resource envelopes in scheduling 
search for problems with multi-capacity resources. This 
includes both direct use as a backtracking and termination 
criterion in a constrained based scheduling algorithm for 
multi -capacity resources and the additional development 
of effective variable and value ordering heuristics based on 
resource envelopes. 

References 

J.C. Beck, A.J. Davenport, E.D. Davis, M.S. Fo> . Beyond Contention: 
Extending Texture-Based Scheduling Heuristics, in Proceedings of AAAI 
1997, Providence, RI, 1997. 

A., Cesta, A. Oddi, S.F. Smith, A Constraint -B a ;ed Method for Resource 
Constrained Project Scheduling with Time Wind ows, CMU RI Technical 
Report, February, 2000. 

H. S.F. Cooper Jr., The Loneliness of the Long-E'uration Astronaut, Air & 

Space/Smithsonian, June/July 1 996, available a 

http://www.ai rspacemag.com/ASM/Mag/Inclex/ 1 996/JJ/Slda.htinj 



